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Abstract Consider the planar linear switched system x(t) = u{t)Ax(t) + (1 — u(t))Bx(t), where A 
and B are two 2x2 real matrices, x £ K 2 , and u(.) : [0, oo[^ {0, 1} is a measurable function. In this 
paper we consider the problem of finding a (coordinate-invariant) necessary and sufficient condition 
on A and B under which the system is asymptotically stable for arbitrary switching functions u(.). 

This problem was solved in previous works under the assumption that both A and B arc diag- 
onalizablc. In this paper we conclude this study, by providing a necessary and sufficient condition 
for asymptotic stability in the case in which A and/or B are not diagonalizable. 

To this purpose we build suitable normal forms for A and B containing coordinate invariant 
parameters. A necessary and sufficient condition is then found without looking for a common 
Lyapunov function but using "worst-trajectory" type arguments. 
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1 Introduction 



By a switched system, we mean a family of continuous-time dynamical systems and a rule that determines 
at each time which dynamical system is responsible of the time evolution. More precisely, let {f u : u £ U} 
(where U is a subset of R m ) be a finite or infinite set of sufficiently regular vector fields on a manifold M, and 
consider the family of dynamical systems: 

x = f u (x), xeM. (1) 

The rule is given by assigning the so-called switching function, i.e., a function u(.) : [0, oof— ► U C K m . Here, we 
consider the situation in which the switching function is not known a priori and represents some phenomenon 
(e.g., a disturbance) that is not possible to control. Therefore, the dynamics defined in (1) also fits into the 
framework of uncertain systems (cf. for instance [5]). In the sequel, we use the notations u £ U to label a 
fixed individual system and u{.) to indicate the switching function. These kind of systems are sometimes called 
"n-modal systems" , "dynamical polysystems" , "input systems" . The term "switched system" is often reserved 
to situations in which the switching function u(.) is piecewise continuous or the set U is finite. For the purpose 
of this paper, we only require u(.) to be a measurable function. For a discussion of various issues related to 
switched systems, we refer the reader to [4, 9, 11, 13]. 

A typical problem for switched systems goes as follows. Assume that, for every fixed u £ U, the dynamical 
system x — f u (x) satisfies a given property (P). Then one can investigate conditions under which property (P) 
still holds for x = f u <t)(x), where u(.) is an arbitrary switching function. 

In [1, 6, 9, 12], the case of linear switched systems was considered: 

x(t) = A u{t) x(t), x£R n , {A u } u&u cR nxn , (2) 

where U C R m is a compact set, u(.) : [0,oo[— > U is a (measurable) switching function, and the map u i— ► A u 
is continuous (so that {A u } ue jj is a compact set of matrices). For these systems, the problem of asymptotic 
stability of the origin, uniformly with respect to switching functions was investigated. 
Let us recall the notions of stability that are used in the following. 

Definition 1 For 6 > let B$ be the unit ball of radius S, centered in the origin. Denote by U the set of 
measurable functions defined on [0,oo[ and taking values on U. Given xq £ l n , we denote by 7 X0 ,«(.)(-) the 
trajectory of (2) based in xq and corresponding to the control u(.). The accessible set from xq, denoted by 
A(xo), is 

A{x ) = U„(. )eW Supp(7 XOi „ ( . ) (.)) • 

We say that the system (2) is 

• unbounded if there exist Xq £ R™ and u(.) £ U such that 7a; 0j u(.) (£) goes to infinity as t — > oo; 

• uniformly stable if for every e > 0, there exists S > such that A(xq) C B e for every xq £ B$; 

• globally uniformly asymptotically stable fGUAS, for short) if it is uniformly stable and globally 
uniformly attractive, i.e., for every 61,62 > 0, there exists T > such that 'YxoM-)^) £ Bs 1 for every 
u(.) £ U and every Xq £ B$ 2 ; 

Remark 1 Under our hypotheses (linearity and compactness) there are many notions of stability equivalent 
to the ones of Definition 1. More precisely since the system is linear, local and global notions of stability are 
equivalent. Moreover, since {A u } ue u is compact, all notions of stability are automatically uniform with respect 
to switching functions (see for instance [3]). Finally, thanks to the linearity, the GUAS property is equivalent 
to the more often quoted property of GUES (global exponential stability, uniform with respect to switching) , 
see for example [2] and references therein. 

Let us recall some results about stability of systems of type (2). 

In [1, 12], it is shown that the structure of the Lie algebra generated by the matrices A u : 

g = {A u : u £ U} l .a., 

is crucial for the stability of the system (2). The main result of [12] is the following: 
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Theorem 1 (Hespanha, Morse, Liberzon) If g is a solvable Lie algebra, then the switched system (2) is 
GUAS. 

In [1] a generalization was given. Let g = r x s be the Levi decomposition of g in its radical (i.e., the 
maximal solvable ideal of g) and a semi-simple sub-algebra, where the symbol x indicates the semidirect sum. 

Theorem 2 (Agrachev, Liberzon) If s is a compact Lie algebra then the switched system (2) is GUAS. 

Theorem 2 contains Theorem 1 as a special case. Anyway the converse of Theorem 2 is not true in general: if 
s is non compact, the system can be stable or unstable. This case was also investigated. In particular, if g has 
dimension at most 4 as Lie algebra, the authors were able to reduce the problem of the asymptotic stability of 
the system (2) to the problem of the asymptotic stability of an auxiliary bidimensional system. We refer the 
reader to [1] for details. For this reason the bidimensional problem assumes particularly interest. In [6] (sec 
also [10]) the single input case was investigated, 

±{t) = u(t)Ax(t) + (1 - u(t))Bx(t), (3) 

where A and B are two 2x2 real matrices with eigenvalues having strictly negative real part (Hurwitz in the 
following), x G R 2 and u(.) : [0, oo[-> {0, 1} is an arbitrary measurable switching function. 

Under the assumption that A and B are both diagonalizable (in real or complex sense) a complete solution 
was found. In the following we refer to this case as to the diagonalizable case. More precisely a necessary and 
sufficient condition for GUAS was given in terms of three coordinate-invariant parameters: one depends on the 
eigenvalues of A, one on the eigenvalues of B and the last contains the interrelation among the two systems 
and it is in 1-1 correspondence with the cross ratio of the four eigenvectors of A and B in the projective line 
CP 1 . A remarkable fact is that if the system x = u(t)Ax + (1 — u(t))Bx has a given stability property, then 
the system x = u{t){r\A)x + (1 — w(t))(r 2 i3)x has the same stability property, for every t\,t% > 0. This is a 
consequence of the fact that the stability properties of the system (3) , depend only on the shape of the integral 
curves of Ax and Bx and not on the way in which they are parameterized. 

The stability conditions for (3) were obtained with a direct method without looking for a common Lyapunov 
function, but analyzing the locus in which the two vector fields are collinear, to build the "worst-trajectory". 
This method was successful also to study a nonlinear generalization of this problem (see [8]). One of the most 
important step to obtain these stability conditions was to find good normal forms for the matrices, containing 
explicitly the coordinate invariant parameters. 

Remark 2 It is interesting to notice that common Lyapunov functions do not seem to be the most efficient 
tool to study the stability of switching systems for arbitrary switchings. In fact, beside providing sufficient 
conditions for GUAS, like those of Theorems 1 and 2, the concept of common Lyapunov function is useful 
when one can prove that, if a Lyapunov function exists, then it is possible to find it in a class of functions 
parameterized by a finite number of parameters. Indeed, once such a class of functions is identified, then in 
order to verify GUAS, one could use numerical algorithms to check (by varying the parameters) whether a 
Lyapunov function exists (in which case the system is GUAS) or not (meaning that the system is not GUAS). 
In the case of stability of switching systems under arbitrary switchings, this is in general a very hard task. 
For instance for systems of type (3) one can prove that the GUAS property is equivalent to the existence of a 
common polynomial Lyapunov function, but the degree of such common polynomial Lyapunov function is not 
uniformly bounded over all the GUAS systems [10]. 

In this paper we provide a necessary and sufficient condition for GUAS in the nongencric cases omitted in 
[6] . In particular we study the stability of the system (3) assuming that at least one of the matrices (say A) is 
not diagonalizable. In the following we refer to this case as to the nondiagonalizable case. We also assume that 
A and B are both Hurwitz and that [A, B] ^ otherwise the problem is trivial. These conditions arc gathered 
in the following assumption (HO), which is often recalled in the following: 

(HO) A and B are two 2x2 real Hurwitz matrices. We assume that A is nondiagonalizable and [A, B] ^ 0. 
We also study the cases in which the system is just uniformly stable. 

A very useful fact is that the stability properties of systems of kind (2) depend only on the convex hull of 
the set {A u } u£ u (see for instance [10]). As a consequence we have 
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Lemma 1 The system (3) with u{.) : [0, oo[— > {0, 1} is GUAS (resp. uniformly stable, resp. unbounded) if and 
only the system (3) with u{.) : [0, oo[— > [0, 1] is. 

In the following we refer to the switched system with u(.) taking values in [0, 1] as to the convexified system. 
Sometimes we will take advantage of studying the convexified system. 

The techniques that we use to get the stability conditions for (3) under (HO) are similar to those of [6] . How- 
ever new difficulties arise. The first is due to the fact that since A is not diagonalizable then the eigenvectors of 
A and B are at most 3 noncoinciding points on CP 1 . As a consequence the cross ratio is not anymore the right 
parameter describing the interrelation among the systems. It is either not defined or completely fixed. For this 
reason new coordinate-invariant parameters should be identified and new normal forms for A and B should be 
constructed. These coordinate invariant parameters are the three real parameters defined in Definition 2 below. 
One (rj) is the (only) eigenvalue of A, the second (p) depends on the eigenvalues of B and the third (k) plays 
the role of the cross ratio of the diagonalizable case. Section 2 is devoted to the computation of the normal forms. 

Once suitable normal forms arc obtained, we look for stability conditions, studying the set Z where the 
two vector fields are linearly dependent, using a technique coming from optimal control (see [7]). This set is 
the set of zeros of the function Q(x) := det(Ax,Bx). Since Q is a quadratic form, we have the following cases 
(depicted in Figure 1): 

A. Z = {0} (i.e., Q is positive or negative definite). In this case one vector field points always on the same 

side of the other and the system is GUAS. This fact can be proved in several way (for instance building a 
common quadratic Lyapunov function) and it is true in much more generality (even for nonlinear systems, 
see [8]). 

B. Z is the union of two noncoinciding straight lines passing through the origin (i.e., Q is sign indefinite). 

Take a point x G Z \ {0}. We say that Z is direct (respectively, inverse) if Ax and Bx have the same 
(respectively, opposite) versus. One can prove that this definition is independent of the choice of x on Z. 
See Proposition 1 below. Then we have the two subcases: 

Bl. Z is inverse. In this case one can prove that there exists u$ g]0, 1[ such that the matrix uqAx + 
(1 — uq)Bx has an eigenvalue with positive real part. In this case the system is unbounded since 
it is possible to build a trajectory of the convexified system going to infinity with constant control. 
(This type of instability is called static instability.) 

B2. Z is direct. In this case one can reduce the problem of the stability of (3) to the problem of the 
stability of a single trajectory called worst-trajectory. Fixed xq £ R 2 \ {0}, the worst-trajectory *f Xo 
is the trajectory of (3), based at xq, and having the following property. At each time t, j Xo {t) forms 
the smallest angle (in absolute value) with the (exiting) radial direction (see Figure 2). Clearly the 
worst-trajectory switches among the two vector fields on the set Z. If it does not rotate around the 
origin (i.e., if it crosses the set Z a finite number of times) then the system is GUAS. This case is 
better described projecting the system on KLP 1 (see Lemma 7 below). On the other side, if it rotates 
around the origin, the system is GUAS if and only if after one turn the distance from the origin is 
decreased, (see Figure 1, Case B2). If after one turn the distance from the origin is increased then 
the system is unbounded (in this case, since there are no trajectories of the convexified system going 
to infinity with constant control, we call this instability dynamic instability). If j Xo is periodic then 
the system is uniformly stable, but not GUAS. 

C. In the degenerate case in which the two straight lines of Z coincide (i.e., when Q is sign semi-definite), one 

see that the system is GUAS (resp. uniformly stable, but not GUAS) if and only if Z is direct (resp. 
inverse). We call these cases respectively C2 and CI. 

The main point of the paper is to translate conditions A., B. and C. in terms of the coordinate invariant 
parameters appearing in the normal forms. Several cases (parametrized by the parameter k) should be studied 
separately: 

• k = 0, called singular case (S-case in the following), 
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Z^ 
(inverse) 

uniformly stable but not GUAS 
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GUAS 



GUAS 



Figure 1: 



• k ^ and B has non real eigenvalues, called "regular —1" case (R_i case), 

• k 7^ and B has real noncoinciding eigenvalues called "regular +1" case (Ri case), 

• k 7^ and B is not diagonalizable, called "regular 0" case (Rq case). 

The case in which B has two coinciding eigenvalues and it is diagonalizable is not considered since in this case 
B is proportional to the identity and, therefore, (HO) is not satisfied. 

The structure of the paper is the following. In Section 2 we compute the normal forms for A and B in which 
the coordinate- invariant parameters appear explicitly. In Section 3 we state the stability conditions, that are 
proved in the next sections. 

Remark 3 As in the diagonalizable case (since the way in which the integral curves of Ax and Bx are 
parametrized is not important), if the system x = u(t)Ax + (1 — u(t))Bx has a given stability property, 
then the system x = u(t)(T\A)x + (1 — u(t))(i~2B)x has the same stability property, for every t\,t% > 0. 

In section 4 we start by studying the set Z where the two vector fields are linearly dependent. In Section 5 we 
state and prove some general stability conditions (in particular cases A, Bl, CI, C2 of Figure 1). In Section 
6 we build the worst-trajectory (i.e., we study case B2). 



2 Basic Definitions and Normal Forms 

For i£R define 

f +1 if x > 
sign(x) = < if x = 

[ -1 if x < 0. 

Definition 2 Assume (HO) and let S be the discriminant of the equation det(B — Xld) = 0. Define the following 



■5 




' if 6^0 



P = 



invariant parameters: 

#^ = o, 

Vj& (5) 
Tr(B) 1 j 

— 5 — V ° = °. 

A ^(AB) - ^Tr(A)Tr(B)) 

7Y(AB) - -Tr(A)7V(B) if 6 = 0. 

Remark 4 Notice that 5 = (Ai — A2) 2 G K, where Ai and A2 are the eigenvalues of B. Notice moreover that B 
has non real eigenvalues if and only if S < 0. Finally observe that 77, p < and k G M. 

Definition 3 7n f/ie following, under the assumption (HO), we call regular case (R-case for short), the case 
in which k ^ and singular case (S-for short), the case in which k = 0. 

Lemma 2 (R-case) Assume (HO) and k 7^ 0. Then it is always possible to find a linear change of coordinates 
and a constant r > such that A/t and B/t (that we still call A and B) have the following form: 

B= ( P sign(S)/k \ ^ 



Moreover in this case [A, B] ^ is automatically satisfied. 

Remark 5 In the following, for the regular case, we call Ri-case, R_i-case, Ro-case the cases corresponding 
respectively to sign(S) = 1, sign(d) = — 1, sign{5) = 0. See Lemma 4 below for the discussion of the eigenvalues 
and eigenvectors of B in these three cases. 

Lemma 3 (S-case) Assume (HO) and k = 0. Then 6 > and it is always possible to find a linear change of 
coordinates and a constant r > such that A/t and B/t (that we still call A and B) have the following form, 

A = ( 7 J ^ I , B = \ ^ ? 1 I 1 called Si -case, (9) 

lo 1/ V p + 1 / 
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A=[ T1 1 ), B=[ P t l 1 ° ), called S-x-case. (10) 



or the form, 

A-f 

7] J ' ^ /O — 1 

Proof of Lemma 2 and Lemma 3. We can always find a system of coordinates such that, 

A =( a)' S= (c d)' a < 6 < c < deR - ( n ) 
In this case the discriminant of B is 5 = (a — d) 2 + 4bc. We have 

m=(; d _7)- (12) 

Notice moreover that c = Tr(AB) - ±rr-(A)Tj-(.B). Hence, according with (6), c = iff k = 0. 
Case c 7^ 0. First notice that in this case [A, B] ^ 0. Consider the transformation 

T=(l f). (13) 



Then 



/ a+d _5_ \ 

A':=r 1 AT = A, £?' := T~ 1 BT = y 2 o|d J • (14) 

we have the following: 

• If <5 = 0, according to Definition 2, then c = k, A = ry, p = (a + d)/2, and with the change of notation 
B' — > i? we get the normal forms (7) and (8). 



If 5 7^ 0, we define 



- A— p= -2- \ „ / a+d_ Vl£[ 

^=4^' = ^ , 5":=4=5'= 71 - . (15) 



1*1 VI- 5 ! 



Consider the transformation 



V2 



I* = ( T ,1 ) • (16) 



We have 



A—-— 1 \ / a-\-d 



4r 



A'" := (T')-^"T' = ( VH ^ , i?'" := [T')~ l B"T' = ^ Q+ ~ d | . (17) 



According to Definition 2, and with the change of notations A 1 " — > A, £?"' — > we get the normal forms 
(7) and (8). Lemma 2 is proved. 

Case c = 0. In this case (HO) implies S = (a — d) 2 > 0. Consider the transformation: 



1 



-6 



T= \ o V )■ (18) 



We have 

4' — T 1 ^ 1 AT 1 — A R' — T^ 1 RT — I 

v d 

We have two cases 



A' := T~ X AT = A, B' := T^BT = ( JJ J ) . (19) 
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• d > a. In this case from a + d = Tr(B) and d — a = \/S, we get 

a= ^(Tr(B)-VS), d = ^(Tr(B) + VS). 

Define 

A" — —A' - I 
A - V6 [ 

Define now T' as in (16). Then 

A'" := (T') _1 A"r 

With the change of notations A 1 " — ► A, — > £> we get the normal forms (9). 

• Similarly if d < a, we get the normal forms (10). 
Lemma 3 is proved. 

The following Lemma can be directly checked. 

Lemma 4 Under the condition (HO), we have the following: 

• A has a unique eigenvalue n < and its corresponding eigenvector is (1,0). 

• if k ^ and S > (i.e., in the Tli-case), then the eigenvalues of B are p + 1 and p — 1, corresponding 
respectively to the eigenvectors (l,k) and (1,—k), and we have p < — 1. 

• if k ^ and S < 0, (i.e., in the R-i-case) then the eigenvalues of B are p + i and p — i, and we have 
p < 0. In this case the integral curves of the vector field Bx are elliptical spirals rotating counter-clockwise 
if k > ( clockwise if k < 0). 

• if k ^ and 5 — (i.e., in the Ho-case), then B has a unique eigenvalue p < 0, and it corresponds to 
the eigenvector (0, 1). 

• ifk = 0, (i.e., in the S-case) then the eigenvalues of B are p+l and p—1, with p < — 1. 

— In the case of the normal form (9) (i.e., in the S\-case) they correspond respectively to the eigen- 
vectors (0,1) and (1,0). 

— In the case of the normal form (10) (i.e., in the S-i-case), they correspond respectively to the 
eigenvectors (1,0) and (0,1). 

3 Main Results 

In this section we state our stability conditions. First we need to define some functions of the invariants ?/, p, 
k defined in Definition 2. Set 

A = k 2 - 4npk + sign(S)4n 2 . 
By direct computation one gets that if k = 2rjp then A = — 4det^4dct_B < 0. It follows 
Lemma 5 Assume (HO). Then A > implies k ^ 2np. 




B" := -^B = 







(20) 



B'" := (T') _1 S"T = B". 



(21) 
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Moreover, when A > and k < 0, define 




l( -k+VA ) 2kp 2 -s l gn(S)(Tr(AB) + VA) l ( VA , g X j£ _ ■ / « 2 / q 

\y- k -VA> 2WdetB(p-sign(,6)%) 1 eXPl k + P°»tgn{6)h 11 P fc ^ U, 



(22) 



exp(- i ^A _|_ pQ_ 1 ), if p — sign(S)j — (which implies sign{5) = — 1). 



where 



arctan 



if fc(p+f) + ?7 ^o 

if fcO>+jJ)+TJ = 0, 



fc(p+f)+»? 



tt/2 



arctanh 



a/A 



fc(p-f)-'?' 



0, 



kp 



Notice that when k < 0, then k(p — ^) — ?y > and fcp > 0. Hence 6\ and are well defined. 

The following Theorem states the main result of the paper. The letters A., B., and C. refer again to the 
cases described in the introduction and in Figure 1. Recall Lemma 5. 

Theorem 3 Assume (HO). We have the following stability conditions for the system (3). 

A. If A < 0, then the system is GUAS. 

B. If A > 0, then: 

Bl. if k > 2rjp, then the system is unbounded, 
B2. if k < 2r]p, then 



C. If A = 0, then: 

CI If k > 2r/p, then the system is uniformly stable (but not GUAS), 
C2 if k < 2ryp, then the system is GUAS. 

4 The set where the two vector fields are parallel 

As explained in the introduction, the stability conditions arc obtained by studying the locus in which the 
two vector fields Ax and Bx are linearly dependent. This set is the set Z of zeros of the quadratic form 
Q(x) :— Det(Ax, Bx). Set x — (xi,X2). We have the following: 



The discriminant of this quadratic form is the quantity A defined in (22). Hence we have the following: 

• if A < then Z = {0}, 

• if A > then Z is a pair of (noncoinciding) straight lines passing through the origin, 





|) x\ + kx\X2 + (r]k)x1 in the R-case 
i + 2nxi) in the Si case 

> — 2r\x\) in the S_i case. 



(23) 
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• if A = then Z is a single straight line passing through the origin (the two straight lines of the previous 
case coincide). 

Remark 6 Notice that, under the condition (HO), in the regular case, the sign of A can be positive or negative, 
while in the singular case, since k = and 8 > (cf. Lemma 3), we have A = 4i] 2 > 0. In the singular case Z 
is the union of a pair of straight lines one of them coinciding with the x\ axis (as it is clear from the expression 
of Q). 

In the following, under the assumption A > 0, we give the explicit expression of the angular coefficients of the 
two (possibly coinciding) straight lines whose union is Z. 

R_x-case (i.e., the case in which B has non real eigenvalues) 



if p + ii/k ^ then m = —. r—r = -—, 7-- , (24) 

if p + i]/k = then m + = 00, m~ = —77. (25) 

Ri-case (i.e., the regular case in which B is diagonalizablc and has real eigenvalues) 



•f /Wn + i ± - k± ^ -k±y/k*-4r,k(p- V /k) 
rfp-77/MOthen m = 2{f> ^ = ^^Jk) ' ^ 

if p — i]/k = then m~ = 00, m + = —Tj. (27) 

Rp-case (i.e., the regular case in which B is nondiagonalizable) 



-k ± VA -fc ± v/fc 2 - 



2p 2p 
S-case (i.e., the singular case) 



(28) 



2?7 

m + = 0, ?77 _ = — < in the Si-case, (29) 

p+1 

2n 

m + = 0, mT = — > in the S_i-case. (30) 

P- 1 

The following Proposition says that if A > and the two vector fields have the same versus on a point of 
Z \ {0}, then this is the case all along Z \ {0}. More precisely it says that if k < 2-qp (resp. k > 2i]p) then they 
have the same (resp. opposite) versus. Recall Lemma 5. 

Proposition 1 Assume (HO), A > , and let Z = D+ U D~ where D ± = {{h,m ± h) el 2 , he E}. Let us 
define a by Bx = a Ax for x € D ± \ {0}. Then 



dot B 



>0, 



dct A 
4- 1 — 2rip—k 

In other words sign(a + ) = sign(a~) = sign(2r\p — k). 

Proof. Let us start with the regular case and A > 0. In this case it is easy to check that a ± = fc ^" . Hence 

a+a- = fc2+fcp(m ^+J+ pm+m ~ ■ Using the fact that m++m" = and m+m" = a5, where X := p-sign(5)%, 
we get 

_t_ _ k(x — p) + p 2r q p 2 rj — sign(S)r] p 2 — sign(S) det B 

oc ct = = = = . 

?7 3 ?7 3 7? 2 det^4 
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Figure 3: 



Similarly we have 

k(m + + m~) + 2pm + m~ —k 2 + 2prjk 2np — k 



a + a = 



nrr&mr rj 2 k det A 

In the regular case, with A = 0, we have D + = D~ , m + = m~ =: m° and a + = a~ =: a . Now Bx\ = oPAx\ 
implies that 

= k + pm° = k + p^ = p -2 X 

Since A = k 2 - Ank X = 0we have 2 X = ^ and a = = Thus 

2a o = 2w~fc and (a o )2 = P 2 signjS) = detS 
det ^4 ?7 2 det 4 

In the singular case k = 0, we have A > and 

-2ry 

Z = {(xi,X2) I x 2 = 0} U {(xi,X2) | x 2 = — — -x\} or 



p+1 
2?] 

Z = {(xi,X2) I x 2 = 0} U {(xi,a52) I ^2 = ——xi\. 

P ^ 



An easy computation show that 



+ _ dctB ,_ 2P 2r 1 p 2 V p 

or a. = 5 = 7 and or + a = — = — — = - -. 

?y z det A V V det A 



Definition 4 If A > and fc < 277/0 (resp. k < 2r\p) we say that Z is "direct" (resp. "inverse"). 

The following Proposition says, roughly speaking, that on Z \ {0} the vector field Ax points "clockwise" (see 
Figure 3). 

Proposition 2 Assume (HO) and A > 0. Let dg be the vector field defined by dg{x\,X2) = (— X2,X\). Then 
in the normal forms of Section 2, we have Ax ■ dg < for every x G Z \ {0}. 

Proof. Let x + G D + , where D + is defined in Proposition 1. Then x + = (h,m + h) for some h G K \ {0}. It 
follows Ax + ■ dg(x + ) = —(m + ) 2 h 2 . Similarly if x~ = (h,m~h) for some h G R\ {0}, we have Ax~ ■ dg(x~) = 
-{m-) 2 h 2 . m 
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5 General Stability Conditions 



In this Section we study some basic stability conditions. 

The following proposition (from which it follows A. of Theorem 3) expresses the idea that if the two vector 
fields are linearly dependent just at the origin, then the system is always GUAS. For the proof one can follow 
the same steps used in the case of two diagonalizablc matrices [6] . It is also a particular case of the same result 
proved for nonlinear systems in [8] . 

Proposition 3 Assume HO. If A < then the system is GUAS. 



5.1 The inverse case 

The next proposition (from which it follows Bl. of Theorem 3) says that if the vector fields are parallel along 
two (noncoinciding) straight lines and point in opposite directions, then there are trajectories going to infinity. 

Proposition 4 Assume (HO) . If A > and k > 2rjp then the system is unbounded. In the degenerate case in 
which A = and k > 2rjp then the system is uniformly stable but not GUAS. 

Proof. To prove the first statement, recall that we can work with the convexified system. Consider the 
following constant control 

oT — \{pi + + a~) 
1 + a + a~ — (a + + a~) 

From Proposition 1 we have a + + a~ < and a + a~ > 0, which implies that uq g]0, 1[. A simple computation 
show that: 

det M («o) = dct( UQJ 4 + (1 - uo)B) = AA , An ,+_ , + T _^ < 0, 

4dct^4(l + a+a — (a + + a ) 

since A > 0. Hence M(uq) has a positive real eigenvalue and the system is unbounded. 

To prove the second statement, observe that if A = and k > 2r]p, then Z = D + = D~ = {(x 1 ,X2) | X2+2rjxi = 
0}, and thanks to Proposition 1, a + = a~ = a < 0. Consider again the matrix 

(a ) 2 - a a 

M(uo) = uqA + (1 — uq)B where uq = 

Then Tr(AI(u )) = ^ozj^ < and detM(u ) = 0. Therefore for u a the system is not asymptotically stable. 
To prove that the system is uniformly stable we show that Ax and Bx admit the following common (non strict) 
Lyapunov function, 

xf 
if] 

Let us prove that 

1 



V{x) = V{x 1 ,x 2 )=x\ + -±. 



VV{x)Ax = ^-(x 2 + 2r ] x 1 ) 2 <0 (31) 
Zr] 

VV(x)Bx = -^{x 2 + 2r / x 1 ) 2 <0. (32) 
2rj z 

To prove (31) observe that W(x)Ax = 2r\x\ + 2x\X2 + = ^j( x 2 + ^V x i) 2 - To prove (32), observe that 
VV { x)Bx = 2 P x\ + + ±_ )xiX2 + J_ xl = 2pA + { ^ 9 nW + P )xiX2 + ^ 

But A = k 2 — Ankp + Ari 2 sign{8) = implies that k 2 + Arj 2 sign{5) = Ankp. Then 

VV(x)Bx = 2px\ + ^£-XiX 2 + = 2 P X 1 + — x l x 2 + 7r^ x 2 = 7To( x 2 + ^"<]X\) 2 . 

2r\ i k 2rf i] 2rf 2rf 

It follows that V(x) is a common Lyapunov function and the system is uniformly stable. 



Notice that in the S-case, Proposition 4 never apply because Z is always "direct" . 



12 



5.2 The direct case: the worst-trajectory 



In the case in which A > and k < 2np, the stability of (3) can be reduced to the study of a single trajectory 
called the "worst-trajectory" . 

Definition 5 Assume (HO) , A > and k < 2r\p. Fix xq £ M 2 \ {0}. The worst-trajectory j Xo is the trajectory 
of (3), based at Xq, and having the following property. At each time t, j Xo {t) forms the smallest angle (in 
absolute value) with the (exiting) radial direction (see Figure 2). 

As explained in the introduction, (see also [6, 10]) the worst-trajectory switches among the two vector fields 
on the set Z. The following Lemma (whose proof is a consequence of the arguments used in [6]) reduce the 
problem of stability of (3) to the problem of the stability of the worst-trajectory. 

Lemma 6 Assume (HO), A > 0, and k < 2np. Fix x G M 2 \ {0}. Then the system (3) is GUAS (resp. 
uniformly stable but not GUAS, resp. unbounded) if and only if "f Xo tends to the origin (resp. is periodic, resp. 
tends to infinity). Moreover when j Xo is periodic or tends to infinity then it rotates around the origin switching 
an infinity number of times. In the degenerate case in which A — 0, and k < 2np, then the system is GUAS. 

Remark 7 The last statement can be be proved either by using worst-trajectory type arguments or by building 
a common quadratic Lyapunov function and it implies statement C2. of Theorem 3. 

Remark 8 The explicit construction of the worst trajectory in the case A > and k < 2np (i.e., the proof 
of case B2. of Theorem 3) is done in the next section. The worst-trajectory can rotate or not around the 
origin. As it will be clear next, if B has complex eigenvalues then it always rotates. One of the statement of 
Lemma 6 is that if the worst-trajectory does not rotate then the system is always GUAS. Notice that, if the 
worst-trajectory rotates around the origin, then it rotates clockwise, thanks to Proposition 2. 

6 Construction of the worst-trajectory 

In this section we prove B2. of Theorem 3, separately for the cases R_i, Ri, Ro, and S. 
6.1 The R !-case 

In the R_i-case we have 5 < 0, therefore the matrix B has non real eigenvalues. 

Notice that under the conditions 8 < and < k < 2r\p we have A = k 2 — 4r/k(p + n/k) = k(k — 2ijp) — 
2knp — Arj 2 < 0. Therefore, since we arc interested to the case in which A > 0, we may assume k < 0. Recall 
formulas (24) and (25). We have the following cases depicted in Figure 4, for x\ > 0. 

• k < and p + n/k > 0. In this case -\/A < |fc|. It follows that m + > mT > 0. 

• k < and p + n/k < 0. In this case ^J~K > |fe|. It follows that m~ > > m + . 

• k < and p + w/k = 0. In this case m + = oo, mT = —n > 0. 

Fix xq = (l,m + ) G D + . In all cases, the worst-trajectory 7^,, rotates clockwise (cf. Remark 8) and is a 
concatenation of integral curves of Ax (from the line D + of equation x 2 = m + x\ to line D~ of equation 
x 2 = m~Xi,) and integral curves of Bx otherwise. See Figure 5. 

Let t\ be the time needed to an integral curve of Ax to steer a point from D + to D~ . Since 




(33) 



we have 



h = 



m 



+ 



m 




13 



R_ t case, A>0, k<2r|p 




k<0, p+r|k>0 k<0, p+riyk<0 k<0, p+ri4c=0 
The worst trajectory rotates around the origin 



Figure 4: 





+ 




yc Bx a 




Ax \( Ax 




\/ Bx ^V 






x m 



Figure 5: 
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Let £2 be the time needed to an integral curve of Bx to steer a point from D to D + . Since 

cost 



y ' vr ' y ksmt cost 

we have 

f arctan fc(p ^f )+? , if fc(p + f ) + r? ^ 

\f if fc(p+£)+r / = 0. 

Notice that in the case p + = 0, we have the following simple formulas for the switching time: f x = — and 
t2 = arctan — . Fix xo G -D + and let 1Z be the ratio between the norm of j Xo after half turn (i.e., after time 
t\ + £2) and the norm of 7 Xo (0), i- c -i 

v _ \\7x (ti + t 2 )\\ 

Il7x„(0)|| [6) 

We have the following cases 
• p + t 7^ 0. In this case 



"ft = cxp(— )exp(pt 2 )|(l H — — )(cost 2 — sint 2 )|, 



after simplification we get 



i.e. 



7 ' =I( ^-Vl ) 2 fc VdetB(p + i) |6XP( — )6XP(pt2) ' 

n = i( fc + ffi 2fcp2 z^lffl^) ± ^ I ex P (^ + 

-/c-VA 2kVdcTB(p- sign(S)l) k 

• p + j: — 0. In this case 

_ /VA X , w m~ . — 2r? /V 7 ^ 

7e= exp — — )exp (pt 2 cosi 2 + — smi 2 ) = = exp(— — + pO-i). 

k k ^Jk 2 + if k 

The system is GUAS, uniformly stable (but not GUAS) or unbounded respectively if TZ < 1, 1Z = 1, 1Z > 1. 
6.2 The Ri-case 

In this section we study the Ri-case i.e., the case in which the matrix B has real eigenvalues. In this case since 
8 > 0, we have A = k 2 — 4r)k(p — rj/k). Again we restrict to the case A > 0, and k < 2r/p. Set 

X--=p-\- (35) 
Recall formulas (26) and (27). We have the following cases depicted in Figure 6, for x\ > 0. 



i) k < 0. In this case we have \ < and yA > \k\. It follows that m~ > > m + and we have: 
Claim 1. k < m + < m~ < —k. 
Proof of Claim 1. We have 

m+ > k <s=^ ~ fc + ^ > k —k + VA < 2\k VA < k(l + 2%). (36) 
2x 
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Ri case, A>0, k<2r|p 




k<0 

the worst-trajectory 
turns around the origin 




k>0, p-T|/k>0 
GUAS 





Figure 6: 



Since 1 + p < (cf. Lemma 4), we have fc(l + 2\) = k(l + 2p — 2-qk) = fc((l + p) + p - 2-qk) > 0. Hence 
m+ > k A < k 2 (l + 2 X f -kr)kx < ixHk + xfy -r] < k + X k < k(p + 1). (37) 
The last inequality holds since p + 1 < 0. Similarly we prove that mT < —k. 

ii) k > and x > 0. In this case \fK > \k\. It follows that m~ < < to + and we have: 

Claim 2. mT < — k < m + < k. 
Proof of Claim 2. We have 

m + > k <^=> ~ fc+ > k ■*=>■ -k + \/A > 2xk ^> V~A > k + 2k\ A > fc 2 (l + 2 X ) 2 

^=> -Ar]kx > 4xk(k + \k) -r) > k + kp - r] ^=^> > k(p + 1). 
Since {p + 1) < 0, the last inequality is always true. In a similar way we can show that to" < —k. 

iii) k > and x < 0. In this case \fK < \k\. It follows that m~ > m + > 0. 

iv) k > and p — rj/k = 0. In this case m~ = oo, to + = — r\ > 0. 

In the case i), the worst-trajectory j Xo rotates clockwise around the origin (cf. Remark 8) and it is a 
concatenation of integral curves of Ax (from the line D + of equation x 2 = rn + xi to line D~ of equation 
x 2 = m~X\) and integral curves of Bx otherwise. See again Figure 5. Let t± be the time needed to an integral 
curve of Ax to steer a point from D + to D~ . Recall formula (33). We have 

to + — m~ vA 

Let t-2 be the time needed to the an integral curve of Bx to steer a point from D~ to D + . Since 

cMBt)=cMpt)[ ^nht c2hi )' 
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we have 

t<i = arctanh 



i.e., 



Defining again 1Z as in formula (34), we have similarly to the previous case 

\J /\ tti^~ yj /\ 

TZ = cxp(— — ) cxp(pt 2 )|(l H ; )(coshi 2 H r~ sinht 2 )|. 

k kr\ k 

After simplification we get 

^'^V^ 2 fc v^(p-i) |CXP( ~ )CXP( ^ } ' 

K = |( ^ + ^) ^ ~ fjWW^) + ^ | cxp( ^ + 
-fc-v^ 2fcVdetS(p-sz 5 n( ( 5)2) k 

As in the previous case, the system is GUAS, uniformly stable (but not GUAS) or unbounded respectively if 

n < 1, ii = i, ii > i. 

Cases ii), iii), iv), are easily studied projecting the system on RP 1 . This is the purpose of the next sec- 
tion. 

6.2.1 The system on the projective space 

In the Ri-case, it may happen that worst trajectory goes to the origin without rotating. In this case the system 
is GUAS. The description of this fact is contained in the following Lemma, where we project the system on 
RP 1 using its linearity. Recall Lemma 6. 

Lemma 7 Assume (HO), A > 0, <5 > 0, k < 2rjp. Consider the projective space RP 1 represented by the 
semicircle {(a;i,2; 2 ) € R 2 : x\ + x\ = 1 and x\ > 0} with the points (0,1) and (0,-1) identified. On RP 1 
represent the eigenvector of A as the point p a = (1,0), the eigenvectors of B as the points 



and the set Z n RP 1 as the couple of point, 



+ I 1 m+ \ - ( ^ m 

Pm ~ ^/T+T^F yiTT^F' Pm VITT^F JT+JmFf 

See Figure 7. Let Ca be the connected component to pa of the set RP 1 \ {p+ ,p~}. If Pb or Pb belongs to Ca, 
then the system is GUAS. 

Sketch of the proof of Lemma 7. By contradiction, if the system is not GUAS, then the worst trajectory 
rotates around the origin and it switches an infinity number of times on the set Z. Under the hypotheses of 
the Lemma, once a trajectory enters the cone 

K A = {xEM. 2 \{0}: ^eC A } (40) 
\x\ 

it cannot leave it. Indeed, if in Ka the worst trajectory corresponds to the vector field B (resp. A), it cannot 
cross the eigenvector of B (resp A) lying in Ca, since pa, P% and p~^ are stable points for the system projected 
on RP 1 ). Hence this trajectory cannot rotate around the origin. Contradiction. H 



In the cases ii), iii), iv), the system is GUAS thanks to Lemma 7. 
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Figure 7: 



B 



R case 
o 




A>0, k<2r|p =^> k<0 



Figure 8: 



6.3 The R Case 

In this section we study the Ro-case, i.e., the case in which the matrix B is also nondiagonalizablc. In this 
case A = k 2 — Ar/kp. Again we assume A > 0, and k < 2ijp. The set Z is the union of a pair of straight 
lines with angular coefficients given by formula (28). Notice that A > 0, and k < 2r\p implies that k < and 
m + < < rri~ . Indeed if k > 0, then A = k 2 — Akrjp = k(k — 2rjp) — 2r\pk < 0. From Lemma 4 we know that 
B has the unique eigenvector (0, 1). The relative position of m + , mT and of the eigenvectors of A and B is 
shown in Figure 8. 

In this case the worst-trajectory rotates clockwise around the origin and it is a concatenation of integral curves 
of Ax from the line D + = {(xi,X2) € R 2 | X2 = m + x\} to the line D~ = {(xi,^) G M 2 | X2 — m~xi}, and 
integral curves of Bx otherwise. In this case we have 



exp(Bt) = exp(pt) 



and the switching times are respectively: 

ti 



1 

kt 1 



— — and t 2 = - — . 

r\k kp 
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A=4r) 2 k=0 



case s case 




A=B= m+ 




A=B= m + 



GUAS GUAS 



Figure 9: 



Defining again 1Z as in formula (34). we have similarly to the previous cases 

.,-k + VE.2kp 2 ~sign(5)(Tr(AB) + v f E). .y/E 

--) », rr— ft/ : — 777^ 1 exp(— — + p6» ). 



-k-y/E 2kVdctB{p- sign{5)l) k 

Again the system is GUAS, uniformly stable (but not GUAS) or unbounded respectively if 1Z < 1, 7Z = 1, 
?e > 1. 



6.4 The S-case 

In the singular case, the worst-trajectory is always tending to the origin (see Figure 9). This is due to the fact 
that the straight line D + = {(h, 0) <E M 2 , h £ R} (belonging to Z) coincide with the eigenvector of A and with 
one eigenvector of B. Using arguments similar to those of Lemma 7, one sees that the worst-trajectory can 
never cross D + and it cannot rotate around the origin. Hence the system is GUAS. One can also check that in 
this case V{x) = x\ + x\ is a common Lyapunov function. 
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